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THE EMISSION-LINE SPECTRA OF MAJOR MERGERS: EVIDENCE FOR SHOCKED OUTFLOWS 
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ABSTRACT 

Using a spectral decomposition technique (Soto & Martin 2012, hereafter Paper I), we investigate the 
physical origin of the high- velocity emission line gas in a sample of 39 gas-rich, ultraluminous infrared 
galaxy (ULIRG) mergers. Regions with shock-like excitation were identified in two kinematically 
distinct regimes, characterized by broad (a > 150 km s _1 ) and narrow linewidths. Here we investigate 
the physical origin of the high-velocity (broad) emission with shock-like line ratios. Considering the 
large amount of extinction in these galaxies, the blueshift of the broad emission suggests an origin on 
the near side of the galaxy and therefore an interpretation as a galactic outflow. The large spatial 
extent of the broad, shocked emission component is generally inconsistent with an origin in the narrow- 
line region of a AGN, so we conclude that energy and momentum supplied by the starburst drives these 
outflows. The new data are used to examine the fraction of the supernova energy radiated by shocks 
and the mass loss rate in the warm-ionized phase of the wind. We show that the shocks produced by 
galactic outflows can be recognized in moderately high-resolution, integrated spectra of these nearby, 
ultraluminous starbursts. The spectral fitting technique introduced in Paper I may therefore be used 
to improve the accuracy of the physical properties measured for high-redshift galaxies from their 
(observed frame) infrared spectra. 

Subject headings: galaxies: starburst — galaxies: evolution — galaxies: active — galaxies: formation 



1. INTRODUCTION 

Ultraluminous infrared galaxies (ULIRGs) are some of 
the most powerful galaxies in the local universe, with 
log(L m /L Q ) > 12. Many ULIRGs exhibit disturbed 
mor phology which indi cates that they are major merg- 
ers (jBorne et al.l [T999). Tidal interactions in major 
mergers are believed to drive gas inflow that fuels both 
the starb ursts and the active galactic nuclei (AGN) in 
ULIRGs (iToomre fc Toomrd Il972t iSanders et al.l 119861: 
ISpringel et al.M2005t Hopkins et al.ll2008D . 

Models of major mergers also predict an early phase 
of sup ernova feedback followed by the emergence of an 
AGN (jSnringel et al.1120051 IHopkins et al.112003 IHookinsI 
120121) . Star burst-driven outflows are found i n 75 — 80% 
of ULIRGs (|Rupke et al.l [20021: lMartin|[2005T) . Interstel- 
lar absorption lines in their spectra are blueshifted a 
few hundred km s _1 relative to th e systemic velocity 
set by the mole cular gas and stars (lArmus et al.l Il990t 
Heckman et al.l 119901: iMartinl I2005L l2006t iRupke et al.l 



2005allbl ld: IRupke fc VeilleuTl2005t ) . AGN-driven out- 
flows, in contrast, have only been found in the small 
subset of ULIRGs with Seyfert 1 nuclei (|Rupke et al.l 
l2005a|) . possibly indicating that these systems are no 
longer classified as ULIRGs by the time AGN drive the 
outflows. Because the measured stellar velocity disper- 
sions of ULIRGs typically exceed the stellar rotation 
speed and th e surf ace brightness profiles fit r 1 / 4 laws 
(|Genzel et al.ll200l . as in elliptical rather than spiral 

1 The data presented herein were obtained at the W.M. Keck Ob- 
servatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California and 
the National Aeronautics and Space Administration. The Obser- 
vatory was made possible by the generous financial support of the 
W.M. Keck Foundation. 



galaxies, removal of the gas - whether by star forma- 
tion, black-hole fueling, or outflow - provides the last 
step in transforming ULIRGs into field elliptical galaxies 
(|Dasvra et al.ll2006h . 

The high SFRs common among z ~ 2 — 3 galax- 
ies are thought to b e fueled by steady gas accretion 
instead of mer gers jDaddi et al.l 120071 : iShapiro et al.l 
2008; INoguchil [L999lV but mergers may play a piv- 
otal role in the f ormation of galactic spheroids even 
at these redshifts (Hopkins et al.l[2011h . Regardless of 
how these galaxies get their gas, comparable star for- 
mation rates are only f ound among ULIRG s and t he 
Lyman-Break Analogs (lOyerzier et al.l 120081 l2009bHal . 
12011 iHeckman et"atl f20lD in the local universe, so 
these environments provide the best local laboratories 
for studying the feedback processes that shape the evo- 
lution of high-redshift galaxies. 

Given the importance of incorporating feedback from 
massive stars and AGN into galaxy formation simula- 
tions, simultaneously mapping out the excitation and 
gas kinematics in these extreme, local environments is of 
broad interest (Heckma n et al.lfl987l: lArmus et al.l Il990t 
Heckman et al.lll990l ; iMurrav et al.l 12005: Vci lleux et al.l 
2009). The emission-line spectrum of shocked regions, 



for example, is easily distinguished from the spectrum 
of gas photoionized by massive stars but quite similar 
to the spectrum of gas photoionized by an AGN. Be- 
cause the size of the narrow-line region ( NLR) powered 
by AG N is limited by the AG N luminosity (|Bennert et al.l 
l2006allbl : iGreene et al.l 1201 lh , one strategy for breaking 
this degneneracy between excitation mechanisms is to 
resolve the location of the gas emitt i ng the shock-like 
spectrum (iMonreal-Ibero et all 120061 : iRich et al.1 120111 : 
IGoncalves et al.1 120101 . Paper I). Large velocity dis- 
persions provide another way to identify regions with 



2 



strong feedback from star formation and active nu- 
clei, gravitational i nstability, and/or streams of recently 
accreted mater i al (iGenzel et al.l 120111 : lLaw et aTl 120091 : 
ISchreiber et al.1 120061 ) and the Doppler shifts of broad 
emission-lines can be used to distinguish shocks gener- 
ated by galactic winds and infall from either cold streams 
or tidally stripped gas. 

To provide a more comprehensive examination of the 
location and speeds of shocks in ULIRGs, we mapped the 
emission- line ratios across 39 local (z = 0.043 — 0.152) 
ULIRGs in velocity and one spatial direction. These 
longslit spectra obtained with the Keck Echellette Spec- 
trograph and Imager (ESI) lack the complete spatial cov- 
erage of integral field spectra (IFU) but offer more sen- 
sitivity and broader spectral coverage than IFU observa- 
tions. In Paper I, we identified the regions with shock- 
like line ratios. Following a brief summary of the sample 
and observations in Section [5] of this paper, the velocity 
dispersion and Doppler shift of the shocked regions are 
presented in Section [3l where we argue that outflows are 
the primary origin of those shocked regions with broad 
(cr > 150 km s _1 ) emission lines. In Section 21 we dis- 
cuss the origin of the outflows and estimate the outflow 
properties. Finally, we discuss the impact of shocked gas 
on the integrated spectrum in Sect. 15.11 

2. OBSERVATIONS AND REDUCTIONS 

This study examines the measurements presented in 
Paper I of 39 ULIRGs at various st ages of merging chose n 
from from the IRAS 2 Jy survey (jMurphv et al.lll996P i. 
These local galaxies span z = 0.043 — 0.163, allowing spa- 
tially resolved optical spectroscopy with the Keck ESI 
l" x 20" long sli t. The data were reduced as described 
in iMartinl (|2005l ). Emission lines in the galaxies have 
~70 km s _1 spectral resolution and ~ 0.8" spatial res- 
olution limited by the atmospheric seeing. The median 
spatial resolution is 1.5 kpc, with a range of 0.7 to 2.1 
kpc depending on redshift. 

3. RESULTS 

In Paper 1, we decomposed the emission line profiles 
into multiple Gaussian components. An example of this 
decomposition is shown in Figure 1. Because we jointly 
fit several transitions, we measured diagnostic line ratios 
for each velocity component. Using kpc-scale apertures, 
we mapped the Doppler shift, velocity width, and line 
fluxes across the ULIRGs. The results, provided in Table 
3 of Paper 1 , were used to classify the excitation mecha- 
nism of each component as Hll-like or shock-like. Here, 
we explore the relationship between the gas kinematics 
and the excitation mechanism. 

For each velocity component at a given spatial lo- 
cation, we classified the excitation mechanism as HII- 
like or shock-like u sing diagnostic ratios of line fluxes 
(jKewlev et al.l 120061) . The [OI] /Ha ratio plays a major 
role in this categorization due to its high sensitivity to 
shocks. In cases where [OI]/Ha is obscured by telluric 
absorption, we rely on [SII]/Ha and [Nil] /Ha to make 
this distinction. If the errors in the flux ratio cross the 
maximum star formation line, then the excitation mech- 
anism is designated "unclear" . We note that our "un- 
cl ear" category d iffers from the composite classification 
of iKewlev et al.l (|2006f) . defined by line ratios between 
the empirically defined maximum star formation limit 
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Fig. 1. — Example of emission-line fit from paper I; the line pro- 
files of the different transtions were fit with common kinematic 
components of varying flux. This figure demonstrates the clear 
variations in the line profiles for each transition, i.e. the blucshiftcd 
wing on Ha, the nearly triangular line profile of [N IIJA6583, and 
the very low intensity blue wing on [O III]A5007. Using this fitting 
method, we are able to deconstruct the line profile and examine 
the excitation mechanism for all of the spectral components that 
show the same kinematics. 

and the maximum star formation. We further describe 
these categories and the fitting method in Paper I. 

Figure 2 shows a histogram of the velocity dispersions 
of the components from all apertures across all galax- 
ies divided by excitation class. Spectral components 
classified as "Hll-like" are narrow in linewidth, having 
a v < 150 km s _1 . The components identified as "shock- 
like", however, present a high- velocity tail at a v > 150 
km s^ 1 . Clearly the broad components are preferen- 
tially shock excited. 

We show the Doppler shift of all the components in Fig- 
ure [3] The narrow, shock-like components (and the HII 
components) are detected over distances of 15 — 30 kpc 
along the slit. The velocity offsets relative to the galaxy 
redshift are positive in one direction and negative to the 
opposite side of the nucleus. Figure 5 of Paper 1 demon- 
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Fig. 2. — In the above histograms we present the distributions 
of line width for components identified as "shock-like" and those 
identified as "HII". Nearly all components with "HII" line ratios 
have a v < 150 km s _1 . The distribution of components identified 
as "shock-like" has a tail out to large line widths. (Of the 430 
components measured, 34 copomnents are not shown here because 
the measurement uncertainties did not place them uniquely in one 
of these two excitation classes.) 



strates that the position- velocity diagrams for these nar- 
row, shock-like components are typically consistent with 
the projection of a galactic rotation curve, suggesting 
the process producing narrow, shock-like emission com- 
ponents is associated with a gas disk. 

The broader components with <j v > 150 km s _1 , how- 
ever, do not show similar rotation gradients, suggesting 
that a different physical mechanism is involved in pro- 
ducing the broad emission. Figure [3] shows that the 
broad, shocked components are nearly always blueshifted 
in stark constrast to the narrow components. The 
blueshifts of the broad, shocked components reach a max- 
imum of 500 km s _1 . 

Because the ULIRGs are very dusty galaxies, their nu- 
clei are not transparent at optical wavelengths. Hence, 
the emission-line radiation that we detect at the cen- 
ter of a ULIRG must be emitted on the near side of 
the dusty, gaseous disk. We therefore conclude that the 
blueshifted emission comes from outflowing gas on the 
near-side of the galaxy rather than infall on the far side. 
The broad emission in adjacent apertures maintain sim- 
ilar kinematics, which allows us to assume that these 
regions have similar dust obscuration. In support of this 
outflow interpretation, we note that the broad, shocked 
emission spans a velocity range quite similar to that of 
the blueshifted Na I a bsorption troughs pr eviously iden- 
tified in these spectra (|Martinll2005l I2006D . 

Using the spatial information along the slit, we traced 
the extent of each spectral component in Paper 1. Figure 
|4] shows how the spatial extent of an emission component 
varies with its velocity disperion and excitation classifi- 
cation. Broad, shocked emission in ULIRGs is usually 
spatially extended but is not detected as far away from 
the nucleus as is the narrow component. The broad, 
shocked component can typically be traced to ~ 5 kpc 
away from the peak continuum emission and up to 6 kpc 
in IRAS 03158+4227. We discuss the source of these 
broad, shocked outflows in Section [4] 

4. OUTFLOW SOURCE 

In general, the broad emission components are 
blueshifted to high velocity relative to the systemic ve- 
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Fig. 3. — We find that the components identified as shock-like 
with a y < 150 km s have a similar distribution in Doppler shift 
as the HII components. The shock-like components with cr v > 
150 km s —1 , however, present significant blueshifts. 
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Fig. 4. — HII and shock-like components with a v < 150 km s 
extend to regions 15 kpc away from their associated nuclei. Com- 
ponents with <r„ > 150 km s _1 , however, remain within 5 kpc of 
the nucleus. The gap between the points at r = r nuc and rest of 
the distribution reflects the spatial resolution of the spectroscopy. 



locity and only a few components are redshifted as shown 
in Fig. [3J For emission, a blueshifted line profile can be 
interpreted in two ways - outflow on the near side of the 
galaxy or inflow on the far side of the galaxy. In the case 
of dusty ULIRGs, however, the presence of dust obscu- 
ration allows a unique identification of the blueshift as 
outflow; photons coming from infalling gas on the far side 
of the galaxy would have a longer path length through 
dusty regions, thereby increasing the amoung of extinc- 
tion. The outflow interpretation is further supported by 
the spatial extent of the blueshifted lines identified in 
Fig. |4j The spatial extent of broad blueshifted emission 
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suggests that the emission is associated with regions be- 
yond the deepest part of the gravitational potential well. 
In the following sections, we explore the possible physical 
origins of the outflow. 

4.1. Star Formation vs. AGN as a Power Source 

The power source of ULIRGs has been under debate for 
more than 20 years, owing to the obscuration in optical 
and shorter wavelengths by dust. Using the Spi tzer In- 
frared Spectrograph (IRS). iVeilleux et all ([2009) exam- 
ined the power source of 14 of the 39 ULIRGs included 
in the Paper I sample by measuring indicators of AGN 
contribution from infrared indicators. Indicators of the 
strength of AGN contributions come from ratios of mid 
infrared emission lines [O I V] / [N e III and [Ne Vl /[Ne I I] 
(lArmus et all l200i I20061 . [20071: IVeilleux et all l2009h . 
however, these ratios are well constrained for only 3 
of the 14 ULIRGs in this sample. Two of the objects 
with well constrained contributions (IRAS01572+0009, 
IRAS12072-0444) show - 90% AGN contributions to 
Lboi- The third object (IRAS15001+1433) that is well 
constrained shows an ~ 50% AGN contribution. Mea- 
surements of the mid-IR line ratios for the other 11 
ULIRGs place a < 30% upper limit on the AGN con- 
tribution to their bolometric luminosity. Although the 
AGN is generally not the dominant source of luminosity 
in these ULIRGs, we still need to investigate whether the 
broad lines could be directly related to the AGN. 

4.1.1. AGN Broad Line Region 

A broad line region (BLR) is a common feature in 
Seyfert 1 and 1.5 galaxies, appearing as permitted and 
intercombination emission lines broadened to 1 , 000 t o 
25,000 kms- 1 FWHM (IQsterbrock fc Ferlandl l2006h . 
The weakness of forbidden transitions in the BLR is ex- 
plained by the high electron density, n e > 10 9 cm" 3 , 
which enhances collisional de-excitation (Peterson 2006). 
Reverberati on mapping places the BLR size at much less 
than 1 pc (jBentz et al.1 120091 : iBrewer et aD 120111 ). sug- 
gesting fast moving gas reflects the deep gravitational 
potential in the vici nity of a supermassive black hole 
([Dietrich et~aLlll999D . In our ULIRG spectra, emission 
from an AGN BLR would be confined to the central aper- 
ture (~ 1"). 

In the ULIRG spectra, however, the broad emission 
lines extend beyond a spatial resolution element, calling 
into question the BLR as a source of these features. Ad- 
ditionally the broad emission features appear not only 
in the permitted transitions, but also in the forbidden 
transitions. These observations, along with the low elec- 
tron density measurements < 100 cm" 3 from [S II] line 
ratios, rule out an AGN BLR as the source of the broad, 
shocked emission components. 

4.1.2. Narrow Line Region 

In the emission-line ratio diagrams, we use the line 
defining the maximum exci tation via extreme star for- 
mation (|Kewlev et al.l 12006) to distinguish HII emission 
from shock-like emission. Above this, line ratios are typi- 
cally attributed to narrow line regions of Seyfert galaxies 
or LINERs. The Seyfert and LINER line ratios can be 
described by fast shocks where the gas density ahead 
of the shock determines the influence of the radiative 



precursor on the overall emission line ratio. Gas rich 
Seyfert galaxies have a larger contribution from the ra- 
diative precursor, resulting in larger [O III]/H/3 ratios. 
LINERs are suspected to have lower gas densities ahead 
of the shock, which leads to lower [O III]/H/3 ratios 
([Dopita & Sutherland! ri995f) . 

The size of shocked AGN NLR region is empirically 
related to the lumin o sity of the [O III] A5007 emission 
([Bennert et al.ll2002t iGreene et all 1201 ID . The scaling 
relation we use (|Greene et al.l 120111 ) is developed with- 
out corrections for in situ dust extinction of Lpm], ow- 
ing to t he difficulty of rob ust corrections for obscured 
quasars (|Reves et al.l [2008). The [O III] luminosities 
(Lr ni]) used in this scaling relation, only includes spec- 
tral components where log([0 III]/H/3) is greater than 
0.5 to avoid flux contributions from an HII region. For 
each galaxy, we estimate Lroini from the broad, shock- 
like components and predict the expected size of an AGN 
NLR. 

Since any nucleus can host an AGN, we measure lumi- 
nosity for each galaxy in the pairs of merging galaxies - 
increasing the total galaxy count in this study to 48. Of 
the galaxies with components that have log([0 III]/H/3) 
> 0.5, the range of luminosities spans several orders of 
magnitude (1.9 x 10 37 < L [0 ni] < 3.2 x 10 42 erg s" 1 ) 
with a median Lroni] = 8 x 10 39 erg s . For 9 galax- 
ies (in 7 ULIRGs - IRAS00188-0586, IRAS01003-2238, 
IRAS05246+0103, IRAS08311-2459, IRAS09583+4714, 
IRAS13451+1232, and IRAS15130-1958) the measured 
outflow radii are within a factor of 1.6 of inferred radii 
from the scaling relation. 

The majority of galaxies have L[om] and spatial distri- 
butions that do not suggest excitation via AGN. In the 8 
ULIRGs (13 galaxies), no components appear above the 
log([0 III]/H/3) > 0.5 cutoff, but we still detect shock- 
like line ratios up to 6 kpc from the closest nucleus, sug- 
gesting some other mechanism is responsible for the ion- 
ization. In 14 of these galaxies the components with 
appropriate log([OIII]/H/3) appear only in regions out- 
side the nuclear aperture. Galaxies that have a strong 
extended AGN presence would more likely consist of a 
contiguous region that includes the nucleus, suggesting 
that these are also not part of an NLR. The remaining 
galaxies have regions with shocked gas at radii greater 
than the estimated NLR size by a factor of 2. 

In general, most cases do not favor AGN as the source 
of the observed emission line profile and instead are con- 
sistent with shocks driven by stellar winds and super- 
novae, as indicated by the AGN fractions. With this 
understanding of the outflow source, we can make esti- 
mates of the outflowing masses and the energy injected 
into the system by these processes. 

4.2. Energetics of Supernova Feedback 

Excluding the AGN-dominantcd cases 
IRAS01572+0009 and IRAS12072-0444, the star 
formation rates in these ULIRGs indicate enormous 
amounts of mechanical power is deposited in the ISM by 
massive stars. In this section, we compare mechanical 
power to the power in the bulk outflow of warm-ionized 
gas. We stress that our estimates of the mass and 
kinetic energy in the bulk flow are no better than factor 
of ~ 3 accuracy. Uncertainty about the filling factor 
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TABLE 1 

Outflow Parameters 
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IRAS00153+5454 


12 


.10 


270 


3.71 


13.31 


51.5 


17 


47 


1.3 


3.6 


68 


0.02 


0. 


.05 


IRAS00188-0856 


12 


.33 


660 


1.57 


2.32 


407. 


130 


31 


58 


13 


128 


0.45 





.10 


IRAS01003-2238 


12 


.25 


1140 


3.66 


3.13 


203. 


67 


155 


21 


50 


96 


0.22 


0. 


.52 


IRAS01298-0744 


12 


.29 


240 


11.27 


44.98 


86.6 


28 


620 


0.6 


14 


57 


0.01 





.24 


IRAS03158+4227 


12 


.55 


1330 


7.22 


5.30 


181. 


59 


423 


11 


80 


191 


0.06 


0. 


12 


IRAS05246+0103 


12 


.05 


250 


2.71 


10.54 


33.1 


10 


22 


1 


2.1 


60 


0.02 


0. 


.04 


IRAS05246+0103 


12 


.05 


850 


3.10 


3.56 


520. 


170 


110 


48 


31 


60 


0.80 





.52 


IRAS08311-2459 


12 


. 10 


600 


4.94 


8.08 


338. 


110 


320 


14 


40 


135 


0.10 





.30 


IRAS09111-1007 


11. 


.98 


290 


1.38 


4.68 


26.3 


8 


7.3 


1.8 


1.6 


51 


0.04 


0. 


.03 


IRAS09583+4714 


11 


.98 


350 


4.42 


12.27 


413. 


140 


320 


11 


26 


51 


0.20 





.50 


IRAS10378+1109 


12 


.23 


1000 


3.80 


3.72 


267. 


88 


150 


24 


11 


92 


0.26 


0. 


.15 


IRAS10494+4424 


12 


.15 


730 


4.34 


5.78 


80.2 


26 


88 


4.6 


15 


76 


0.06 





.20 


IRAS10565+2448 


11 


.98 


250 


2.66 


10.54 


23.2 


8 


29 


0.7 


2.8 


66 


0.01 





.01 


IRAS11095-0238 


12 


.20 


260 


3.65 


13.61 


145. 


18 


120 


3.5 


8.8 


95 


0.04 





.09 


IRAS12071-0444 


12 


.31 


570 


7.95 


13.57 


15.9 


5 


210 


0.1 


15 


25 


0.02 





.61 


IRAS15130-1958 


12 


.03 


890 


3.65 


4.03 


29.5 


10 


51 


2,1 


13 


24 


0.10 





.53 


IRAS16090-0139 


12 


,18 


640 


6.61 


10.07 


17.0 


6 


86 


0.6 


8.5 


179 


0.003 


0. 


.05 


IRAS16487+5447 


12 


.12 


400 


5.79 


14.01 


147. 


19 


360 


3.5 


26 


71 


0.05 





.36 


IRAS17028+5817 


12 


.11 


200 


1.33 


6.42 


121. 


10 


11 


6.2 


2.1 


69 


0.09 


0. 


.03 


IRAS18368+3549 


12 


.19 


370 


4.09 


10.70 


333. 


110 


110 


10. 


13 


83 


0.12 


0. 


.16 


IRAS19297-0406 


12 


.36 


730 


4.30 


5.74 


233. 


77 


180 


13 


31 


123 


0.11 





.25 


IRAS20087-0308 


12 


.39 


330 


5.62 


16.48 


1170 


390 


760 


21 


16 


132 


0.18 


0. 


.35 


IRAS23327+2913 


12 


.03 


630 


3.81 


5.95 


60.5 


20 


99 


3.4 


17 


58 


0.06 


0. 


.29 


IRAS23365+3604 


12 


.13 


580 


0.89 


1.50 


5.52 


2 


1.6 


1.2 


1.1 


73 


0.02 





.02 



Note. — Column 1: IRAS Name. Column 2: Total infrared luminosity logfLiR/L©) flMurphv et al.iri996f) . Column 3: Outflow velocity taken 
as u ou t — |f| + <7 V ( km s — 1 ) from the Paper I (Table 3). Since we are unsure of the degree to which projection effects play a role in the measured 
velocity, assume all outflowing regions have the same u ut, but is foreshortened. In this simplification we take the maximum i> ut of the apertures, 
to represent the vector that is most parallel to the line of sight. Column 4: Maximum radius of outflow (kpc). r ou t is defined by the maximum 
separation between an aperture which shows outflow and the continuum peak. Column 5: outflow time scale (Myr). Column 6: Total luminosity 
of Ha from the outfl ow components (Lt o t/10 40 ) erg s _1 . We correct the Ha emission for Gal actic dust extinction using the reddening curve from 
ICardell i ct al. (1989) and values of E(B -V) and Ay from the Infrared Science Archive maps fSchlcgcl ct al. 199J§). We use the measured Balmcr 
decrement and the lCalzetti et all d 19941) reddening curve to correct for internal extinction. Column 7: Mass of outflowing gas in 10 6 Mq from Ltot 
and Eq[^ assuming n e — 100 cm -3 and 7' — jphot- Column 8: Mass of outflowing gas in 10 6 M from £ and Eq|6]with an assumed hemisphere 
geometry, with the maximum path length through the sphere defined by the radius, Ss — r out , and / — 0.01. Within each aperture the path length 
through the sphere varies as s — (r^ ut — (^aperture — Nucleus) 2 ) 1 '' 2 Column 9: Outflow rate estimated from Ltot in yr _1 described in Section 
14.2.21 Column 10: Outflow rate estimated from S in Mq yr — 1 with an assumed hemisphere geometry. Column 11: Star formation rate estimated 
from Lir are values JMurphy ct al. 1996) and relation SFR — Lift/13.0 X 10 9 Lq. This relation is an adjustment of the Kcnnicutt ( 1998) relation 
adjusted for a stel lar mass range from 1 to 100 Mq (Mq yr~ 1 ). This star formation rate is also corrected for AGN contribution to L D oi for the 
objects included in Vcillcux ct al. (2009). In cases where upper limits to the AGN fraction are not provided for an individual object in [Veilicux et all 
(2009). wc use the average AGN contribution of 30%. Column 12: Efficiency of outflow from Ml. Column 13: Efficiency of outflow from Ms with 
an assumed hemisphere geometry and filling factor / — 0.01. The objects in this table are restricted to those that have at least 1 aperture with a 
shock- like spectral component with a v > 150 km s~ . 



of the warm-ionized gas and the large-scale geometry 
of the outflow dominate the systematic error. The 
mass, energy, and momentum carried by the outflowing, 
warm-ionized gas should clearly be viewed as a lower 
limit on the total amounts carried by all phases of the 
wind, which likely includes substantial components of 
both coronal gas and molecular gas. 

4.2.1. Mass of Warm- Ionized Gas 

For each ULIRG that shows broad, shock-like emis- 
sion components, we can estimate the outflowing mass 
of warm-ionized gas. This mass is the product of the gas 
density and the volume occupied by the warm-ionized 
phase. 

The fraction of the total volume V filled with warm- 
ionized gas is called the filling factor. For HII regions, 
this filling factor, 



is computed from the RMS electron density obtained 
from the Ha luminosity or surface brightness and the 



electron density computed from a density-sensitive dou- 
blet ratio such as [S II] AA6717, 31. Values of the filling 
factor range from 0.001 < / < 0-1 in nearby H II re- 
gions l|Searleil97lUKennicudll984HKaufman et al.lll987| ) 
and therefore have a very signficant impact on the esti- 
mated mass. The range of filling factor may indeed be 
different in the ULIRGs, given the difference in star for- 
mation rate from these HII regions in normal galaxies, 
but these values provide a point of comparison. In the 
ULIRG spectra, the flux ratio for the [S II] doublet is typ- 
ically consistent with the low-density limit and indicates 
n e < 100 cm~ 3 , but the ratios of the two [S II] lines indi- 
cate densities n e « 500 cm~ 3 in the central apertures 
of IRAS01003-2238, IRAS05246+0I03, IRAS083f f-2459, 
IRAS09583+47I4, and IRAS23327+29I3. For our mass 
estimates, we adopt a fiducial value of n e ~ 100 cm~ 3 
but carry along the scaling with density explicitly in our 
results. 
The outflow mass is 

M out = fimiiYn e f, (2) 

where n e is the electron density, inn is the mass of hy- 



() 



drogen, and the mass per H atom is /i = 1.4 amu when 
helium is included. The unknown factor Vn e f in Eqn. [2] 
is simply related to the Ha luminosity by 



(3) 



Eliminating the common factor between Eqn. 2 and 3, 
the outflowing mass, 



/xm H L Hc 



(4) 



can be estimated from the luminosity of the broad, 
shocked emission component for any assumed value of 
the gas density, n e . 

The effective volume emissivity, 7', varies depend- 
ing on the temperature and excitation mechanism. For 
photoionized regions with T = 10 4 K, 7' = 7 P hot = 
a^hz; = 3.56 x 10~ 25 erg cm 3 s -1 in Cas e B recombi- 
nation theory (Os terbrock fc Ferlandl [2006h . In shocked 
regions, collisional excitation and ionization can make 7' 
larger or smaller than 7 D hot by a facto r of two to three 
(jGoerdt et all 120101: IGenzel et alll20Tl . Since the pho- 
tionization case is in the middle of the plausible range 
for 7 P hot, we estimate masses using 7p hot- 

We present the estimated masses in Table [TJ They 
range from 1.8 x 10 6 M Q to 3.9 x 10 s M . The median 
and standard deviations are 40 x 10 6 M Q and 84 x 10 6 M© 
respectively. The systematic errors in the mass estimate 
are difficult to quantify. For example, Fig. [5] illustrates 
the aperture correction to the observed flux. In addition, 
corrections for internal Ha extinction can be large. The 
large outflow mass in IRAS20087-0308, for example, is 
due in large part to the extinction correction. For indi- 
vidual galaxies, the masses in Col. 7 of Table 1 should 
be interpreted as rough estimates accurate to a factor of 
2-3. 

The measured Ha surface brightness provides a sanity 
check on the assumed electron density, or equivalently 

(rig) 1 ^ 2 Z^ 1 / 2 by Eqn.[T] The Ha surface brightness, £ = 
Lho./A depends on the three-dimensional shape of the 
outflow since the volume is the product of the projected 
area A and the depth of the emitting region, As, along 
the sightline. The surface brightness depends on density 
as 

£ = (n 2 ) 7 'As (5) 



By analogy to Eqn. 4, the inferred outflow mass is 



Ah 



( limn AT, 



As/ 

7'S 



(6) 



In Table 1, col. 8 provides examples of the mass estimates 
obtained by this method. 

If our density estimate n e is reasonable, then equat- 
ing Ms and Mr, should imply reasonable values for the 
product of the filling factor and the depth of the emitting 
region along the line of sight. If we adopt / = 0.01, then 
the required path length through the emitting region is 
often only a few hundred pc. Such a small pathlength 
would arise only if the emitting region had a hollow 
shell structure. If the emitting region is roughly spher- 
ical, As « y/A, the implied filling factor is often small, 
/ ~ 10~ 3 . The real values of s and / that unify the mass 




Fig. 5. — The long black rectangle represents the slit, the gray 
boxes represent apertures along the slit, the annuli are the total es- 
timated regions for each position. The number of o utflows detected 
in ULIRGs (75-80% IHiipke et al.ll2u03 ; IMartinll2005l) suggests that 
a biconical outflow cone has a large opening angle ~145°. This 
allows us to use a sphere as a rough approximation to the geome- 
try. We estimate the luminosity contributed by regions outside of 
the slit by assuming this circular symmetry and a constant surface 
brightness within the half annuli as the aperture within it. The 
central aperture coincides with the nucleus of the galaxy used in 
the measurement. 

estimates may involve a change in both of these param- 
eters, but will depend on each individual case. The bot- 
tom line is the volume filling factor of the warm-ionized 
gas may be tiny thereby leaving lots of volume to be filled 
by gas at a different temperature. 

4.2.2. Outflow Velocity and Time Scale 

The time scale for the outflow is estimated from the 
geometry and the kinematics measured from the emission 
line profile. In Table [IJ we define the outflow velocity 
(w ut) as |i>| + o~ v ( km s" 1 ) taken from Table 3 in Paper 
I. The median v out in this sample is 580 km s _1 , and 
the values range from 200 to 1330 km s" 1 . 

Estimating the outflow timescale from the outflow ra- 
dius and velocity, t out — r out / v OM t , the implied mass loss 

rate of warm-ionized gas is M = M out /r. For the es- 
timate of mass outflow using luminosity (Ml), we find 
a range from 0.4 to 58 M yr -1 where the median is 
4.6 M yr" 1 . The largest M is in IRAS00188-0856, due 
mostly to the moderate v out and small r out . The sur- 
face brightness mass estimate with a spherical geometry 
(Ms) gives larger outflow rates on average from 1 to 80 
M yr" 1 where the median is 15.2 M yr -1 . In this case 
the largest outflow rate is in IRAS03158+4227, which has 
a large v out and r out , but in the surface brightness case, 
the large r out has a greater influence. 

4.2.3. M vs SFR 

An important parameter in modeling galactic winds is 
the efficiency (77) that star forming regions are able to 
remove gas from the star forming region. This parame- 
ter informs the timescale over which a galaxy will form 
stars before the stellar population itself removes gas and 
shuts down star formation. We estimate 77 as the ratio of 
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outflow mass to t he star formation rate measured from 
IR luminosities in iMurphv et al.l (|1996[ ) (Table [J). We 
use the star form ation rate estim ated from the SFR = 
L m /13.0x 10 9 L Q (jKennicuttll 19981 ). adjusted for a stellar 
mass range from 1 to 100 M Q and correcting the Ltr for 
the AGN fraction. 

In this sample of galaxies, the outflow rates typically 
yield r\ ~ 10 _1 but span the range 10~ 3 < 77 < 1. We 
emphasize, however, that this measurement refers only to 
the warm ionized gas, which may be the dominant phase 
of the outflowing gas. The mass outflow rates in the cool 
gas traced by Na I are at least a few tenths of the star- 
formation rate, 77 = M/SFR « 0.1, but are poorly con- 
strained due to the la rge variations in ionization param - 
eter among ULIRGs ([Martini 12005 iMurrav et alil2007[ ). 
Hot winds and molecular outflows may also carry signif- 
icant mass. 

4.2.4. Shock Energetics 

The presence of outflows implies that mechanical en- 
ergy is injected into the gas. For objects that show out- 
flow in emission, we can estimate the energy injection 
rate into the gas from a starburst by using the star for- 
mation rate. We sc ale the predicted sup ernova energy 
injection rates from (Leith erer et al.l 1999, SB99) to the 
calculated SFR for each source, after correcting for differ- 
ences in the assumed IMF low mass cutoff (Kennicutt 0.1 
- 100 M Q , SB99 1-100 M ). We assume continuous star 
formation model over a timescale longer than 40 Myr. 
This produces a nearly continuous injection of feedback 
energy contributed by supernovae and stellar winds. 

The simplest model for the initial response of the inter- 
stella r gas to this feedba c k is based on stellar w ind bub- 
bles (j Weaver et al.lll977t IShuTl fc McKeel fl979T ) . Shocks 
direct 55% of the feedback energy into the thermal en- 
ergy of low density, coronal bubble and the other 45% 
into an expanding shell of swept-up interstellar medium. 
About 60% of the shell energy is radiated away in this 
model, suggesting that luminosity of the shocks would be 
27% of the feedback power. For typical shock velocities, 
most of this energy will come out in ultraviolet and op- 
tical emission lines. The energy in the optical lines will 
be slightly greater than the UV lines by ~ 20%, suggest- 
ing that ~ 15% of the feedback power will come out in 
optical emission (IShull fc McKedll979D . 

Using our estimates of physical properties in Section 
14.2.11 and the shock speeds from Paper I, we can make 
comparisons to these efficiencies. In Figure[S]we show the 
ratio of the total optical luminosity (L opt ) for the compo- 
nents indicating shocked outflowing gas in the measured 
lines compared to the mechanical power (L w , inferred 
from the SFR) . Since the luminosity of the optical emis- 
sion lines is generally less than ~15%, it follows that su- 
pernova feedback is a plausible source of power for this 
shocked emission. 

Luminosity in the optical emission lines for 4 
of the galaxies (IRAS05246+0103, IRAS09583+4714, 
IRAS13451+1232, and IRAS20087-0308) exceeds the ex- 
pected fraction of energy injected into outflow by super- 
novae. The excess emission line energy in these objects 
(except for IRAS20087-0308), may be due to additional 
emission from an AGN NLR, since for these objects, the 
emission region size from scaling relations is consistent 
with that of an AGN NLR. Emission from IRAS20087- 



0308 encounters heavy extinction making the Balmer 
decrement difficult to measure, resulting in a possible 
over correction to the extinction, which would then skew 
the L opt /L w ratio. 

We calculate the ratio of power that appears as kinetic 
energy from the outflow rate (Lke) to the mechanical 
energy injection rate from the estimated mass and veloc- 
ities measured for the outflow regions in Table [T] This 
kinetic energy estimate uses Ml as the outflow mass esti- 
mate and v ou t for the velocity. In Figure [HI we show that 
this energy rate is again consistent with the fraction of 
energy in outflow models (mean = 0.08), which predict a 
~ 20% fraction of mechanical energy in the kinetic energy 
of the shell. The consistency in kinetic energy fraction 
with the models suggests that this process is possibly at 
work driving the outflows in ULIRGs. It also implies 
that the estimated masses are correct within an order of 
magnitude. 
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Fig. 6. — above: A histogram showing the measured luminosity 
(Lopt) emitted in shocked gas compared to the mechanical energy 
injection rate (L w ) from supernovae scaled from the IR star for- 
mation rate with SB99. For a whole galaxy, both the emission 
associated with the shocked outflow (hashed) and the emission 
from all shocked regions (solid) are generally in the range - 15%, 
with the peak below 5%. The measured values are consistent with 
IWeaver et al.l J 19771) ■ below: The fraction of the feedback energy 
that contributes the kinetic energy in the mass outflow is again 
consistent with the Weaver model. 



5. DISCUSSION 



5.1. Outflows in Integrated Spectra 

In observations of galaxies at high redshift, obtaining 
similar spatial resolution to this study is difficult due to 
the small angular size of galaxies. The degree to which 
spatially integrated emission lines show evidence for the 
outflows in these spectra is unclear. We examine spa- 
tially integrated spectra in these data to examine which 
of the resolved features remain evident in the integrated 
spectra. 

5.1.1. Spectral Classifications from the Total Line Flux 

In low spectral resolution, spatially unresolved studies 
of galaxies, one means of examining the excitation mech- 
anism is through total line flux classifications. We apply 
the same classifications (Paper I) to the ULIRGs based 
on spectra obtained by integrating along the slit. These 
total line flux classifications may be influenced by the ob- 
served shocks in outflows and gas disks. We examine the 
relation between classification and the presence of these 
shocks to determine which of the Seyferts, LINERS, and 
HII classified galaxies host shocked gas. 

All 6 galaxies classified as Seyferts present broad, 
blueshifted outflows in the kpc scale apertures. The 
narrow emission line components in these Seyferts also 
present shock-like ratios in four out of six cases - 
IRAS00188-0856, IRAS01003-2238, IRAS13451+1232, 
and IRAS15130-1958. These four Seyfert galaxies 
have both narrow and broad emission components that 
present shock-like flux ratios. The luminosity of the 
broad, shocked component is sufficient Lrom] to generate 
a spatially resovled NLR. We conclude AGN photoion- 
ization is the source of excitation for these four galaxies. 

IRAS09583+4714 and IRAS12071-0444 present the 
two exceptions for the Seyfert classified galaxies, where 
the broad components differ from the narrow compo- 
nents in their excitation. Of the two galaxies in the dou- 
ble IRAS09583+4714, the galaxy with shocked outflow 
presents narrow component excitation consistent with 
HII and has a strong, disk-like rotation gradient. In 
IRAS12071-0444, the narrow components have excita- 
tion more consistent with a LINER classification and are 
highly extended to ±9 kpc. For both of these excep- 
tional galaxies, scaling relations from L[om] suggest that 
the outflow region is larger than what would be expected 
from an AGN NLR. The inconsistency of the outflow ex- 
citation with the integrated classification in these two 
galaxies suggest that a significant fraction of ULIRGs 
may be misclassified as Seyfert galaxies, when they host 
shocks in outflows outflows driven by star formation. 

The subset of 9 galaxies classified as LINERs fre- 
quently host both narrow, shocked components and 
broad, shocked components. In all but one of these ob- 
jects, narrow, shocked emission appears least 5 kpc from 
their associated nucleus. Among the 6 LINER spectra 
with broad, shocked components, the outflow is spatially 
extended in 5 of them. None of the 9 galaxies have suf- 
ficient L[oni] to suggest that the shocked outflow or the 
narrow, shocked component is powered by an AGN NLR. 
These examples emphasize that a integrated galaxy spec- 
trum with a LINER classification need not indicate the 
presence of an AGN. 

Among the sample, 12 galaxies had emission line ra- 
tios on the borderline between classifications. Six of these 



borderline galaxies have ratios between the Seyfert and 
LINER classifications. These six galaxies have broad out- 
flows as well as extended shocks that are not part of the 
bulk outflow. The two borderline galaxies with classi- 
fications between HII and LINER do not host shocked 
outflows, but do have extended shocks, which is respon- 
sible for making the net classification slightly outside the 
range of the HII galaxies, the last of these twelve galax- 
ies are a mix of classifications, without a clear trend in 
behavior. 

The largest fraction of the sample (21 out of 48 galax- 
ies) is classified as HII in the integrated spectra. These 
galaxies share the same HII classification in the diagnos- 
tic diagrams that compare [O I] /Ha and [S II] /Ha, but 
the [N II] /Ha can be ei t her H II or "composite" as de- 
scribed in Ke wlev et al.l (|2006D . The subset of galaxies 
with a "composite" classification are also host to either 
shocked outflows or shocks in the extended narrow emis- 
sion line gas. Clearly, the presence of shocks in these 
objects modifies the [N II] /Ha classification, but obser- 
vations with kpc scale spatial resolution showed that su- 
pernovae are a plausible source of the forbidden line en- 
hancement and the outflow. 

The spatially resolved analysis was also aided by the 
simultaneous multiple component fit to the emission line 
profiles, so we next investigate the integrated spectrum 
with the multiple component fit. 

5.1.2. Identifying Shocks in the Integrated Spectral Profiles 

Our analysis of spatially and spectrally resolved line 
emission provides new insight on how to identify regions 
of shocked gas. Here, we examine whether the method 
works on integrated spectra with moderate spectral res- 
olution. This analysis differs from the previous section 
by comparing the line shapes in the integrated profiles 
to those in the apertures. 

This analysis allows us to define the required ranges 
and resolutions of instruments obtaining the integrated 
spectra. Clearly, the integrated spectra must cover 
[O I]A6300. The strong [O I]A6300 emission from shocked 
regions and the separation of this transition from other 
lines make it critical for uniquely fitting multiple velocity 
components. Furthermore, we will assume the integrated 
spectra have spectral resolution of at least R > 5000 and 
signal to noise ~ 15 if multi-component fitting will be 
attempted. 

Of the 24 galaxies that exhibit broad blueshifted out- 
flow in the spatially resolved analysis, 19 galaxies also 
show these features in the integrated spectrum. In the 
remaining 5 galaxies, the strong narrow lines observed in 
the resolved analysis mask the signature of broad outflow 
in the integrated spectrum. 

Regions with narrow (a v < 150 km s _1 ) emission 
lines that have shock-like line ratios appear frequently 
(34 of 48 galaxies) in the spatially resolved analysis. 
Over a large range in position, these narrow features are 
identified at Doppler shifts ranging from -250 km s _1 to 
250 km s in a few galaxies. Averaging over a large 
range in Doppler shift acts to broaden the integrated 
line profiles. In 5 of the 34 galaxies with narrow shock 
emission, the summed components form a line of width 
<7 V > 150 km s _1 in the integrated spectrum. These 
galaxies are mostly classified as HII with flux ratios in the 
"composite" region of the [N II] /Ha vs [O III]/H/3 diag- 
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Fig. 7. — Plot of the Ha+[N II] and [O I] regions of the spectrum for both the integrated apertures and the spatially resolved apertures 
for IRAS16090-0139 and IRAS23365+3604. Line color represents the excitation classification of each component into "HII" (red), "Shock- 
like" (blue), "unclear" (green). The line style identified the components that are part of an identified outflow (solid), compared to the 
rest of the emission components (dashed). Left: IRAS16090-0139 presents outflow components in the resolved analysis that also appears 
in the integrated spectrum. The observed narrow shock also appears in the integrated profile fit. ffight: In the nuclear aperture of 
IRAS23365+3604 we identify outflow from the kinematics in the emission component. The apertures at 2.9 kpc and 3 kpc from the 
continuum source identify emission from shocked gas. All of these features are washed out in the integrated profile. 
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nostic diagram. Fig. 0shows the integrated aperture and 
spatially resolved apertures for IRAS16090-0139, which 
exhibits the outflow and narrow shocks from the resolved 
analysis in the integrated profile. We also compare this 
to the case of IRAS23365+3604, which shows both types 
of shock in the spatially resolved analysis, while the in- 
tegrated case is only identified as HII. 

In 15 out of the 34 galaxies, the Doppler shifts of the 
narrow, shocked components span a small enough range 
to be detected as 'narrow, shocked' components in the 
integrated spectrum. A narrow span of Doppler shifts 
would also make any broad outflow more identifiable, 
which means that objects with low inclination gas disks 
will show broad outflows more often. 

We conclude that applying the multiple component fit- 
ting technique to integrated spectra would have some 
success in identifying broad, shocked components from 
outflows. In our sample, the outflow is undetected in 
the integrated spectrum only 20% of the time. In these 
5 failures, a broad component is in fact detected in the 
integrated spectrum; but the spatially resolved analysis 
indicates the emission comes from a gas rotating like a 
disk. Building a better understanding of the physical 
origin of the narrow, shock-like emission component is 
therefore critical for interpreting the presence of broad, 
shocked emission in integrated spectra. 

5.2. Lessons From the Local Laboratories 

While there are important differences between ULIRGs 
and star- forming disks at z > 2, ULIRGs present high 
specific star formation rates that are similar to the 
high end of s tarforming galaxies from 2. 5 < z < 5 
(lElbaz et all [20071 ; iCunha et al] l2010t iFeulner et ail 
l2005h . The similar specific star formation rates suggest 
that local ULIRGs can provide insight into processes that 
occur below current resolution limits and that the broad 
emission components may share a common origin. Us- 
ing adaptive optics (AO) spectroscopy of z ~ 2 galaxies, 
star-forming clumps were recently sh own to be the origi n 
of the broad Ha/[N II] line wings (jGenzel et al.M201lD . 
The ~ 4 kpc offset of these clumps from the galactic 
center suggests a supernova-driven outflow, rather than 
an AGN, broadens the line emission. The more central 
location of the broad emission found in the ULIRGs il- 
lustrates that the ULIRGs are not direct analogs of the 
z r~j 2 galaxies but the radius and luminosity of the broad 
emission region in most of these ULIRGs is uncomfort- 
ably large to be caused by the AGN. 

In the echellete spectra, broad emission features are 
identified in the line profiles through the simultaneous 
fitting of multiple transitions. Our fitted linewidths are 
not as large as the broad Hq/JN I I] wings reported for 
z ~ 2 galaxies (Sh apiro et al.l[2009| ) because our spectra 
require broad components in the forbidden lines as well 
as the recombination lines. Additionally, since ULIRGs 
are known to be much more reddened than the z ~ 2 
starbursts; broad emission- line profiles produced by sim- 
ilar physical processes in the less-dusty z ~ 2 galaxies 
will be more symmetric with respect to the systemic 
veloci ty and therefo r e bro ader than those in ULIRGs. 
While (Shapiro et al.l (|2009f) had difficulty distinguishing 
a BLRs and superwinds as a source of the broad emis- 
sion features identified in stacks of z ~ 2 spectra, the 
simultaneous fitting of more broad forbidden lines allows 



the identification of other sources of the broad emission 
rather than just Ha. 

The estimates of outflow mass and outflow rate in the 
local ULIRGs provides context to the estimates of out- 
flow rate and mass in the starforming clumps at z ~ 2 
(jGenzel et al.| [2011). The star formation rate in these lo- 
cal ULIRGs (24 - 180 M yr" 1 ) is much greater than the 
star formation rate in the z ~ 2 starforming clumps (3.3 
to 40 M Q yr _1 ), but our estimates of the mass outflow 
rate are much lower (0.4 to 58 M yr _1 for ULIRGs, 6 
to 200 M Q yr -1 for starforming clumps). These results 
are consistent with a higher rj for warm ionized gas at 
high redshift but do not demand it. The gas density or 
the filling factor adopted for the high redshift galaxies 
could be tuned to give results for 77 similar to what we 
infer for these ULIRGs. 

Some properties of these outflows are easier to study in 
nearby galaxies, and the methods in this study illustrates 
how more information about the outflows at z ~ 2 can 
be obtained, even in integrated spectra. As mentioned 
already, for example, the outflows in z ~ 2 galaxies likely 
have shock-like line ratios similar to those in our ULIRG 
sample. The properties of these shocks can be estimated 
once spectral coverage of additional lines - particularly 
[O I], H/3, and [O III] - are obtained at moderate spectral 
resolution, R « 5000, with similar signal to noise (~ 15) 
as in this study. 

6. SUMMARY 

Paper I mapped optical, emission-line ratios across 39 
ultraluminous infrared galaxies in 2 dimensions, i.e., line- 
of-sight velocity plus one spatial dimension. In this pa- 
per, we have used those measurements to describe the 
kinematics of regions excited by different physical pro- 
cesses. Figure 2 shows the distribution of line-of-sight 
velocity dispersion for regions with shock-like and HII- 
like line ratios. The median linewidth for the shock-like 
components, a v ~ 144 km s _1 , is higher than that of 
the Hll-like components, 61 km s _1 . The difference is 
caused by the large number of shock-like components 
broader than a v f=a 150 km s _1 and the near absence of 
Hll-like components with linewidths this large. 

These broad emission components are relevant to our 
understanding of gas inflows and outflows because gas 
in virial equilibrium would produce lines with smaller 
widths. Our results provide insight into the interpreta- 
tion of the emission-line spectra of high redshift galaxies, 
particularly dusty galaxies with extremely high star for- 
mation rates. 

We find that the broad shock-like components are typ- 
ically spatially extended, reaching radii up to 6 kpc in 
Figure 4. We show that in 9 out of 24 broad components, 
the size of this region exceeds that expected for an AGN 
NLR under the assumption that [OIII] 5007 luminosity 
from components with log([0 III]/H/?) > 0.5 is powered 
by the AGN. We estimate the total power radiated by 
summing the optical luminosities of the broad line com- 
ponents and applying a geometrical correction for slit 
losses. Since the resulting shock luminosities range from 
5-20% of the mechanical power from supernova explo- 
sions, we conclude that supernovae are a viable power 
source for the gas flows that generate the shocks. 

In ultraluminous infrared galaxies, dust absorbs much 
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of the starlight; and the centers of these galaxies are not 
transparent at ultraviolet and optical wavelengths. At 
least in the central few kpc, we can be reasonably cer- 
tain that the emergent emission-line profile is shaped by 
gas on the near side of the galaxy. We therefore interpret 
the blueshifts of the broad components as direct evidence 
that the gas is outflowing. In less dusty galaxies, this out- 
flow component would presumably also have a red wing 
from emission coming from the far side of the galaxy. 

The broad components have surprisingly smooth line 
profiles. Across luminous infrared galaxies with winds, 
for example, the Ha + [Nil] profiles are well described 
by a single, Gaussian component where emission from 
HII regions in the underlying disk dominate the flux. 
When the outflows are observed against the sky, along 
the min or axis, the emission lines are typically double 
peaked (jLehnert fc Heckmanl 119951 ). The relative inten- 
sities of the two components have been shown to reflect 
the inclination and op ening angle of a biconical outflow 
(|Heckman et al.l I199O0 . That some of the broad emis- 
sion is detected beyond the continuum emission in these 
ULIRGs, yet does not have a double-peaked profile may 
be related in part to the obscuration of the back side of 
the outflow. However, it is not obvious that the galaxies 
are opaque at these large radii, and we tenatively con- 
clude that the warm-ionized, outflowing gas in ULIRGs 
does not share the biconical structure that describes out- 
flows in LIRGs very well. 

The prominence of shock-excited emission in these 
starburst galaxies raises questions about how shocked 
emission skews the emission-line ratios measured from in- 
tegrated spectra. Since the luminosity in the ultraviolet, 



ionizing continuum exceeds the mechanical power from 
superno va and stellar winds by about an or der of mag- 
nitude (jLeitherer et al.l 119991 : Martini l2007fl . we would 
not expect shocks from galactic winds to determine the 
line ratios measured in integrated spectra. What Fig- 
ure 8 clearly demonstrates, however, is that when inte- 
grated spectra have moderately high spectral resolution, 
R ~ 5000, shocked, outflowing gas can sometimes be rec- 
ognized by comparing the line profiles of forbidden lines, 
particularly [01] 6300, to the blended Ha + [Nil] profile. 
We conclude that multi-component line fitting of at least 
these 4 transitions in the rest-frame optical spectrum can 
provide useful diagnostics of shock velocities and radia- 
tive energy losses in galaxies over a very broad range of 
cosmic time. 
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